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ABSTRACT 


It is often a challenge to achieve uniform flow in turbulent swirl flow and to predict the flow within the nozzle as measurement 
diagnostics face difficulty to capture both mean flow and turbulence. The purpose of this study is to numerically investigate the near 
wall flow characteristics and turbulent behavior for the effect of different tangential inlet numbers of an incompressible turbulent 
swirl air jet. In this regard, axial-plus-tangential flow based swirling nozzle is considered for the simulation using finite volume 
method, where turbulence is approximated by the Shear Stress Transport (SST) k-co model. The results show that axial and tangential 
velocity at the wall vicinity response the most. Moreover, the turbulent flow characteristic for no swirl flow is nearly uniform, but 
for swirl flow it fluctuates abruptly near the inlet section where the swirl has introduced. The skin friction coefficient for 2TP is the 
maximum for swirl flow and for no swirl condition the skin friction coefficient is nearly uniform. Due to the swirl introduction the 
pressure drop characteristics near the nozzle center response quickly and near the wall vicinity this property changes slowly. The 
magnitude of swirl decay fluctuates before the nozzle converging section however after the nozzle converging section the swirl 
decay is nearly constant. The local swirl near the inlet is highly unpredictable although after the nozzle converging section the local 
swirl profile is nearly similar for 2TP, 3TP and 4TP. 
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1. Introduction 

Swirl flows are a classical form of fluid flow and typically 
encountered in many engineering applications, such as gas 
turbine combustor and cyclonic separator, as well as in nature, 
such as tomedo. In swirling flows generally two types of vortex 
emerge, namely solid body rotation and free vortex flow. 
However, in actual practice, a mixed of these two typically 
appears whereby azimuthal velocity increases with radius in the 
core and decreases afterwards. The academic and industrial 
research on swirling jets was found to be challenging due to their 
unique characteristics, such as flow reversal, vortex bubble and 
anisotropic turbulences in both free and impinging jets. The 
detailed flow physics, instability and above features of swirling 
flows are explicitly available in the literature [1]. The axial and 
radial pressure gradient in swirling flows affects downstream 
flow development and wall-bounded flows. Swirl flows are 
generated in different applications in many ways, such as rotating 
pipe, rotating vanes and twisted tapes inside stationary pipe, 
axial-plus-tangential entry [2]-[7]. The wide variations of 
generations can be divided into two broad groups: geometrical 
and aerodynamical. In geometrically generated swirl flows, for 
example rotating vanes or twisted tapes inside pipes, the flow is 
disturbed by such geometries or obstacles, which exacerbates the 
generality of the flow and difficult to draw a summative 
observation. On the other hand, aerodynamically generated swirl 
flows, such as rotating pipe or axial-plus-tangential entry flow, 
have better control of the flow and common flow physics may 
establish in different researches. Despite geometrically generated 
swirl flows are widely available in the literature and still 
progressing, fundamentals of aerodynamic swirl flows are 
currently being investigated in detail because of its prevalence to 
use recently in number heat treatment applications. 


Yajnik and Subbaiah [8] experimentally studied the effects 
of swirl on internal turbulent flows by conducting experiments 
on turbulent pipe flow. They found that the mean velocity 
profiles close to and away from the wall admit similarity 
representations at sufficiently large Reynolds numbers and 
extended velocity defect law is sensitive to swirl as the wall law 
is not sensibly dependent on swirl. Chang and Dhir [9] found a 
flow reversal region in axial velocity profile in the central portion 
of an acrylic tube and an axial velocity increase near the wall. 
Kito and Kato [10] studied the near wall velocity distribution of 
turbulent swirling flow in circular pipes and concluded that the 
flow becomes three-dimensional after transitional swirl intensity. 
Later, Kitoh [11] investigated the turbulent behavior of free- 
vortex-type swirling flow through a long straight circular pipe 
and reported that the swirling component decays downstream as 
a result of wall friction. Buschmann et al. [12] explored the wall 
behavior, the location of the peak Reynolds shear stresses and the 
three normal stresses of turbulent channel/pipe flows, and stated 
that no scaling works equally for all parameters. Ahmed et al. 
[13] conducted measurements of mean velocity and turbulence 
of swirling flows using dual wire CTA probe both in the core and 
near-wall regions. However, the measurements were confined to 
the immediately above the nozzle exit plane. Effect of viscosity 
and surface tension of fluids and associated instabilities in 
annular flow were also analyzed in some studies [14]-[15]. 

Since both DNS and LES typically require sufficiently 
smaller sized mesh or large mesh quantity to resolve small-scale 
turbulence, Reynolds Averaged Navier-Stokes (RANS) were 
widely used in swirling flows by many investigators. Although it 
is believed that among various RANS models the Reynolds 
Stress Model (RSM) would perform better because of its ability 
to capture anisotropy of turbulent shear stress, this hypothesis 
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was not found to be universal in the past [16]-[17]. In fact, some 
studies revealed that two-equation models (based on Boussinesq 
approximation) can be applied with acceptable accuracy in 
comparison with experimental data for moderate swirl flows - 
[18]-[19]. Nouri-Borujerdi and Kebriaee [20] simulated the 
turbulent boundary layer of an incompressible viscous swirling 
flow through a conical chamber using finite volume method and 
reported that the boundary layer thickness is dependent on the 
velocity ratios, Reynolds number and nozzle angle. Najafi et al. 
[21] investigated turbulent swirling decay in a vertical straight 
fixed pipe where swirl in induced by rotating honeycombs. They 
showed RSM with two-layer zone model for different near wall 
approaches are fairly well to predict the swirling flow but fail to 
predict the pressure distribution along the pipe wall. Islam et al. 
simulated an aerodynamically generated swirl flow using SST k- 
co model in the nozzle exit plane. The result showed centerline 
velocity [22] decay for introducing low levels of swirl into the 
impinging jet and a significant reduction in turbulent kinetic 
energy at the wall region. More recently, researchers [23]-[27] 
performed simulations of aerodynamic swirl flow using SST k- 
cd model and validated their results with the literature. 



(a) 3D Model (b)Diff event Tangential 

Ports 


The above discussion revealed that although a significant 
amount of research available on swirling flows, they are mostly 
confined to gas turbine combustors, cyclonic separator, or 
geometry induced swirl flows. Relevant aerodynamically 
swirling flows are either classical in nature (experimental) with 
outdated experimental facility to resolve near-wall and 
turbulence behaviors inside the nozzle or tested robustness of 
various numerical schemes against experimental data with 
limited focus on flow behaviors. Thus, a detailed study for the 
effect of geometric parameters on near-wall and turbulence 
behaviors for an aerodynamic swirl nozzle appears to be limited 
in the literature. As such, the current numerical study will bridge 
this gap by investigating non-swirling and swirling flows from 
an aerodynamic nozzle for the same initial and flow conditions. 
The paper will examine the effect of number of tangential ports 
on mean and turbulence flow development along the length of 
the nozzle. 

2. Methodology 

2.1 Problem formulation 

An aerodynamically generated swirl nozzle, which is 
capable of seamless transition from non-swirling to highly 
swirling jets, is considered in this study. The nozzle is axial- 
plus-tangential entry type and consists of three tangential 
around the nozzle periphery and an axial port at the bottom of 
the nozzle. Detailed dimension of this nozzle is available in 
[4],[28], hence is not repeated here for brevity. The uniform 
flow was ensured by the flow settling chamber with 
honeycomb. The aerodynamic swirl flow was generated when 
both the uniform axial and tangential flows from different 
number of circumferentially oriented and inclined ports mix 
together farther downstream. In this case, three variants of 
tangential ports: 2-Tangential ports (2TP), 3-Tangential ports 
(3TP) and 4-Tangential ports (4TP) are considered. The exit of 
the nozzle has diameter (. D ) 40 mm, with a total length of 577 
mm. A three-dimensional view with the relative orientation of 
tangential ports is shown in Fig. 1. 


Fig. 1 (a) CAD view of the aerodynamic swirl nozzle (sliced 
to show the internal cavity), and (b): different number of 
tangential ports to impart azimuthal component. 


2.2 Governing Equation and Boundary Conditions 


The governing equations to solve the incompressible, 
steady-state flow characteristics within the aerodynamic swirl 
nozzle are the conservation of mass and momentum as follows: 


v-y = 0 


a) 


pV-VV = -VP + /iVT 


( 2 ) 


Here, ^ is the velocity vector, and p and p are the density 
and dynamic viscosity of the fluid, respectively. Since the 
current problem is turbulent in nature, RANS approach is 
applied to solve the mean and turbulence quantities. In RANS 
approximation, each variable is composed of time-averaged 
part (steady) and turbulence part, as shown below: 


$-(/)+(/)' ^ 

where ^ is a variable used in equations (1) and (2). Upon 
implementing the equation (3) into the governing equations and 
using equation (1) and setting time-average of turbulence 
equals to zero, the resulting RANS equations emerge. The 
RANS equations are similar to the governing equations, except 
an additional term ( pu[Uj ) in equation (2). This additional term 
is known as turbulent shear stress, which governs the 
turbulence characteristics. The shear stress components are 
determined via mean velocity gradients by the Boussinesq 
hypothesis: 


u i u 'j =-kSij-ju t 


' dui_ + du£ 

dx j dx; 

\ J 1 J 


(4) 
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Where, jut is the turbulent viscosity and is a function of k and co, 
which are determined via an appropriate turbulence model. In 
this research, the SST k-co model, proposed by Menter and 
Egorov [29], is chosen to model turbulence transport quantities: 
turbulence kinetic energy, k , and the specific dissipation rate, 
co. SST k -co model was found to be one of the best performing 
RANS models in swirling and wall-bounded flows [16], [30]- 
[31]. 

The boundary conditions considered in this problem are 
mass flow inlet at the axial and tangential ports, pressure outlet 
at the nozzle exit plane and wall at the nozzle wall. The inlet 
conditions are adopted from the study [4] for realistic 
predictions. For no-swirl flow, axial flows (only) are provided 
from the bottom of the nozzle, with no flows from tangential 
ports. Conversely, for the swirling flows, three tangential ports 
supply the same amount of mass flow, with no axial flow from 
the bottom of the nozzle. This ensures the same Reynolds 
number, where the average velocity in Reynolds number is 
determined by volume flow rate divided by nozzle exit area. 
The pressure outlet condition is applied at the nozzle exit with 
atmospheric pressure whereby turbulence is specified by 2% 
intensity and hydraulic diameter. Finally, no-slip condition at 
the nozzle wall with ambient temperature is used. 

2.3 Numerical Settings and Model Validation 

The above governing equations (1-4) and two transport 
equations (SST k-co) are solved by finite volume based 
commercial software ANSYS Fluent vl7. The pressure-based 
coupled algorithm is used to simultaneously solve the 
governing equations. For the pressure discretization PRESTO 
(PREssure STaggering Option) is applied, and the second-order 
central-differenced for diffusion terms and second-order 
upwind scheme for convective terms of the transport equations 
are used. Hexahedral mapped mesh type was used in multi¬ 
zone meshing with an element size 0.0025m. A typical mesh 
independence test is shown in Fig. 2 for 3TP. In this case, four 
different grid elements are tested, namely, 974001, 1310326, 
1647089 and 1892870 elements. The results are found to be 
invariant, except the 974001 elements. As such, a mesh 
containing 1647k element is chosen in this study. Another 
testing is also done for first layer height on y + values and wall 
shear stress profiles, and is shown in Fig. 3. It appears the first 
layer cell height has little influence on those wall 
characteristics. As such, a first layer height of 0.05 mm is 
chosen in this study. The figure also shows that above settings 
ensured y + values less than 0.2 for the whole domain. The 
solution is assumed to be converged when the residuals of the 
flow parameters are less than 10" 5 . Mass conservation is also 
checked between inlet ports and the nozzle exit for converged 
solutions and the difference between inlets and outlet is found 
close to zero. Cartesian coordinates are given by (v, y, z) with 
corresponding velocity components ( U x , U y , U z ). For 
convenience, results of this 3D simulation are presented by 
polar coordinates (v, r, 0) with x-axis coinciding for each 
system. In this regard, the corresponding axial (U) and 
tangential (W) velocity components are derived from the 
Cartesian coordinate data using the axes transformation rule. 



Fig. 2 Sensitivity of mesh elements for the case 3TP. 



(a) 



(b) 


Fig. 3 Effect of first layer cell height on y+ and normalized 
wall shear stress. 


sum of flowrates through all axial and tangential inlets. 
Mathematically, 

Q = Q 

Qa + Qt Or (5) 

The Reynolds number is defined as, 


2.4 Parameter Definition 


Av 


( 6 ) 


The flow ratio Q r is the ratio of total mass flow rate 
through tangential inlets to the total flow rates in the nozzle, i.e. 
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Where, A is the nozzle area at the exit plane of the nozzle 
and v is the dynamic viscosity of fluid. Two different swirl 
number definitions are common as the ratio of axial flux of 
tangential momentum to the axial flux of axial momentum: 


s = - 


2 jcp^ r 2 UWdr 


2 J*Rp\q r\u 2 -0.5W 2 \dr 


JO 


f r~UWdr 

O' — l0_ 

rR 2 

In rU dr 


R 


Finally, the local swirl number is defined as, 


(7) 


( 8 ) 



* 


s 



2.5 Validation 


(9) 


The simulation data of the swirl nozzle is first tested with 
different turbulence models and validated by comparing with 
experimental data [4] at the nozzle outlet plane with three 
tangential ports for the case Q r = 0 and 1 (Fig. 4). It appears 
that SST k-co model predicts the flow behavior well for both 
non-swirling and swirling cases. As such, SST k-co model will 
be used for the data presented in the ensuing results and 
discussion section. Fig. 4 also shows a good agreement 
between the numerical prediction and the experimental data for 
Q r = 0. A slight deviation is observed for Q r = 1, but 
importantly, the numerical data predicts well the profile 
behaviors and peak locations. The deviation is attributed to the 
possible measurement inaccuracy associated with CTA X-wire 
and experimental flow settings [13]. This results an 
overestimation of the mean velocity components than their 
corresponding true values. The deviations may also be partly 
attributed to the inability of RANS approach to accurately 
capture highly swirling flows. 



(a) 



(c) 

Fig. 4 Different turbulence models against experimental data 
[4] for (a) U/U_b for Q_r= 0, (b) U/U_b for Q_r= 1 and (c) 
W/U_b for Q_r= 1. 

3. Results and Discussion 

This section includes presentation of mean and turbulence 
controlling parameters, such as velocity, boundary layer 
thickness, pressure drop, wall shear stress and turbulent shear 
stress for the effect of different number of tangential ports at 
two flow conditions (Q r = 0 and Q r = 1). 

Swirl decay along the length is shown in Fig. 5 for the 3TP 
case. The swirl number is calculated using the Equation (7) and 
Equation (8). It is observed that near the nozzle inlet the swirl 
number S is highly fluctuating and after the converging section, 
the swirl number is nearly constant. Initially, the swirl number 
decreases from the inlet. When the tangential port has 
introduced the value of swirl number increases and 
immediately after the increment the swirl number drops again. 
In the nozzle converging section the magnitude of the swirl 
number rises again at first and then the value decreases and 
finally increases before coming at a constant magnitude. In case 
of S' the values initially increasing from the nozzle inlet come 
to a constant magnitude after the nozzle converging section. It 
is evident that the swirl number S shows a very unpredictable 
nature before the nozzle converging section. 
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Fig. 5 Swirl decay along the length of the nozzle for 3TP case. 

Fig. 6 shows the local swirl numbers along the radial 
location of the nozzle at the different axial positions ( x/D = 
1.275, 3.1, 5.65,9.4, and 14.425) of the nozzle for three 
different numbers of tangential ports, namely, 2TP, 3TP, and 
4TP. It is observed that the profile of local swirl at x/D = 
5.65,9.4, and 14.425 are almost the same in nature for all the 
cases. The local swirl increases gradually from the center of the 
nozzle toward the wall and near the wall, a sudden drop occurs. 
However, the local swirl profiles at x/D = 3.1 are similar for 
2TP and 3TP although for 4TP the local swirl profile is 
parabolic from the nozzle center toward the nozzle wall. The 
value of local swirl near the inlet at x/D = 1.275 is highly 
unpredictable due to the swirl induction. The local swirl profile 
for 2TP and 3TP are nearly similar but for 4TP the profile is 
different. The local swirl near the nozzle wall suddenly 
increases and drops immediately after the increase and finally 
raise to the wall for both 2TP and 3TP. However, the local swirl 
profile at 4TP is like a bell shape; increasing from the nozzle 
center it becomes constant at r/D = 0.2 to r/D = 0.35 and 
then decreases again towards the nozzle wall. As the swirl is 
introduced at this location the local swirl profile becomes 
highly unpredictable. 

The axial velocity vectors at five different axial locations 
(xD = 4.75, 7.50, 9.75, 12.5 and 14.425) are presented in Fig. 
7 for non-swirling (Fig. 7a) and swirling (Fig. 7b-d) 
conditions. Velocity is found to be the almost uniform after the 
converging section for no swirl flow (Q r = 0), but the 
magnitude is the highest near the exit plane, with a reduction 
towards the wall due to boundary layer formation. In contrast, 
for swirl flow ( Q r = 1), velocity magnitude is found to be the 



higher near the inlet and towards the nozzle exit, the velocity 
magnitude decreases. The axial velocity magnitude is 
proportional to the number of tangential ports and the velocity 
magnitude is zero at the wall due to no-slip condition. The axial 
velocity vectors at x/D = 4.75 is symmetric from the 
centerline when swirl is induced (Fig. 7b-d), however, the 
velocity distribution becomes asymmetrical as the flow 
approaches the outlet. The velocity distribution is more 
symmetrical along the radial direction for 2TP than 3TP & 4TP 
with a slight decrease near the center except for exit plane. At 
x/D = 14.425, velocity reduction near the center is highest for 
2TP. 



Fig. 6 Local swirl charactersitics along the length of the 
nozzle (a) 2TP, (b) 3TP, (c) 4TP. 
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Fig. 7 Axial velocity vectors at different axial locations for non-swirling (a), and swirling flows: (b) 2TP, (c) 3TP and (d) 4TP. 
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Fig. 8 depicts the boundary layer thicknesses along the 
length of the nozzle for both non-swirlign and swirling (2TP, 
3TP and 4TP) conditions. It is observed that the magnitude of 
boundary layer thickness at 2TP is pretty low compare to the 
other. For swirl flow conditions the boundary layer thickness at 
the beginning is small and near the exit the boundary layer 
thickness increases. At the middle of the nozzle the boundary 
layer thickness is nearly constant. However, for no swirl flow 
condition the boundary layer thickness higher and then the 
boundary layer thickness gradually decreases towards nozzle 
exit. 



3 6 9 xW 12 15 


Fig. 8 Boundary layer thickness along the length of the 
nozzle. 

Fig. 9 demonstrates the axial ( U/Ub ) and tangential ( W/Ub ) 
velocity components near the wall i.e. at y/D = 0.01,0.05 
and 0.1 for Q r = 0 and Q r = 1. In the figure, Y represents the 
distance from the wall of the swirl nozzle. It is observed that at 
for swirl condition the near-wall axial velocity profile is 
identical for all the location and it increases from the inlet to 
the nozzle converging section. Then, the velocity suddenly 
decreases at Y = 0.0ID immediately after the converging 
section. The axial velocity is then nearly constant for all the 
axial location of the nozzle. However, the axial velocity 
decreases very little after the nozzle convergence at Y = 0.05 D 
and 0.1D, and no significant change occurs. For swirl flow 
condition the near-wall axial and tangential velocity increases 
after the swirl flow introduction. Then after the nozzle 
converging section the near-wall axial and tangential velocity 
show a very little deviation and along the nozzle axial location 
the velocity profile is nearly uniform at Y = 0.05 D and Y = 
0.1 D. However, the near-wall axial and tangential velocity 
profile at Y = 0.0 ID is fluctuates highly along the length of the 
nozzle, especially the tangential velocity profile shows a totally 
unpredictable nature. The overall near-wall axial and tangential 
velocity distribution at Y = 0.0ID is low in magnitude than at 
Y = 0.05D and Y = 0.1 D. It can be observed that the velocity 
profile at the wall vicinity (Y = 0.0ID) response highly than 
the other position. 

The interplay between the number of tangential ports and 
the swirl intensity for the effect of skin friction coefficient 
along the length of the nozzle is shown in Fig. 10. It is observed 
that for no swirl condition the magnitude of Skin friction 
coefficient is very low and along the nozzle axial location no 
change occurs. However, the Skin friction coefficient at high 
swirl flow at 2TP is the maximum along the nozzle axial 
location although for 3TP and 4TP no significant deviation is 


observed. In every case of swirl flow the Skin friction 
coefficient at the nozzle exit is the maximum. Although the 
Skin friction coefficient along the axial location is not 
significantly changing however, near the exit it suddenly 
increase. The overall distribution of the Skin friction 
coefficient is not similar for all the case. 

Fig. 11 displays the pressure drop at different radial 
locations for swirling flow (3TP). The inlet pressure is taken at 
the reference pressure and pressure drop is calculated based on 
the pressure of the nozzle inlet. It is observed that the pressure 



Fig. 9 Near-wall normalized axial (U/Ub) and tangential 
(W/Ub) velocity distrution along the nozzle length: Q r = 0 (a), 
and Qr— 1 (b and c). 

drop at the center of the nozzle (r/D = 0)is greater than all 
other radial location in almost every axial position of the 
nozzle, however, near the nozzle outlet it suddenly reduces. 
While pressure drop at every radial location increase after the 
introduction of the tangential ports, at (r/D = 0.4) it suddenly 
decrease. Pressure drop near the wall is gradually increasing 
after the converging section of nozzle but pressure drop near 
the nozzle center is nearly constant. It can be concluded that 
Due to the tangential flow introduction the pressure drop 
characteristics near the nozzle center response quickly and near 
the wall vicinity this property change slowly. 
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Fig. 10 Skin friction coefficient along the nozzle length. 



Fig. 11 Pressure drop along the nozzle length for swirling 
flow. 



U. 1 U. ~ Urn D T (/. J Q 01 02 t 

r/D r/D 

Fig. 12 Turbulent normal stress profiles at different axial locations. 


Fig. 12 presents the radial distribution of normalized 
turbulent normal stress components at various axial locations. 
It is observed that the turbulent flow characteristic near the wall 
is very high at x/D = 9.4 and x/D = 14.425 at no swirl flow 
( Q r = 0) although at the center of the nozzle the turbulent flow 
characteristic is very low. Moreover, near the inlet of the nozzle 


at x/D = 1.275 and x/D = 4.45 the overall turbulent flow 
characteristic is very small in magnitude. On the contrary at 
swirl flow condition ( Q r = 1) near the inlet of the nozzle at 
x/D = 1.275 the overall turbulent flow characteristic is very 
high. However, the turbulent normal stress near the outlet 
section at x/D = 9.4 and x/D = 14.425 is small in these 
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locations. While at no swirl condition the turbulent normal introduction of swirl flow. The overall magnitude of the 

stress characteristic at the center of nozzle is very low at swirl turbulent normal stress is very small in no swirl flow condition 

flow condition, the u'u' and w'w' component is very high at but for swirl condition it is very high, 

the nozzle center especially at x/D = 1.275 due to the 
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Fig. 13 Turbulent shear stress profiles at different axial locations. 


Finally, Fig. 13 presents the radial distribution of 
normalized turbulent shear stress components at various axial 
locations. It is observed that for no swirl flow ( Q r = 0) the 
turbulent shear stress component u'v' and l/w'near the wall is 
very high at x/D = 9.4 and x/D = 14.425 although at the 
center of the nozzle the turbulent share stress component is very 
low (Fig. 13a and Fig. 13c). However, the u'w' component at 
x/D = 4.45,9.4 and 14.425 is reducing near the wall of the 
nozzle. Near the inlet of the nozzle at x/D = 1.275 the overall 
turbulent share stress is uniform and does not deviate much for 
no swirl flow ( Q r = 0). On the other hand at swirl flow 
condition ( Q r = 1) near the inlet of the nozzle at x/D = 1.275 
all the three component of turbulent share stress is fluctuates 
highly and the magnitude of overall turbulent share stress is 


very high. However, the overall turbulent share stress at x/D = 
9.4 and x/D = 14.425 is very uniform along the radial 
location for swirl flow. The overall magnitude of the turbulent 
flow characteristic is very small in for swirl flow compared to 
the swirl condition. 

4. Conclusion 

An incompressible turbulent swirling air jet is investigated 
numerically. In this regard, an axial-plus-tangential swirl flow 
is considered for non-swirling (Q r = 0) and highly swirling 
(Q r = 1) cases. Finite volume based commercial software 
ANSYS Fluent vl7 is used in the simulation to investigate 
mean flow and turbulence characteristics. Governing equations 
are approximated by the RANS equations and turbulence is 
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characterized by the SST k-co model. The study examines the 
effect of number of tangential inlets on mean flow behaviors 
and turbulent characteristics. The magnitude of swirl decay is 
fluctuates highly before the nozzle converging section and after 
the nozzle converging section the swirl decay is nearly 
constant. The boundary layer thickness in swirl flow is found 
to be the smallest for 2TP. The skin friction coefficient along 
the nozzle axial position at no swirl flow is uniform and the 
magnitude is very small. For swirl flow, the skin friction 
coefficient along the nozzle axial position is fluctuating 
especially for 2TP. Moreover, the skin friction coefficient at the 
nozzle exit is the maximum for all the case. Due to the 
tangential flow introduction, the pressure drop near the nozzle 
center response quickly and near the wall vicinity this property 
changes slowly. The turbulent normal and shear stress for no 
swirl flow is nearly uniform, but for swirl flow it fluctuates the 
most near the section where the swirl flow is introduced. 
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